Conductance of a spin-1 quantum dot: the two-stage Kondo effect 
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We discuss the physics of a of a spin-1 quantum dot, coupled to two metallic leads and develop a 
simple model for the temperature dependence of its conductance. Such quantum dots are described 
by a two-channel Kondo model with asymmetric coupling constants and the spin screening of the 
dot by the leads is expected to proceed via a two-stage process. When the Kondo temperatures 
of each channel are widely separated, on cooling, the dot passes through a broad cross-over regime 
dominated by underscreened Kondo physics. A singular, or non-fermi liquid correction to the 
conductance develops in this regime. At the lowest temperatures, destructive interference between 
resonant scattering in both channels leads to the eventual suppression of the conductance of the 
dot. We develop a model to describe the growth, and ultimate suppression of the conductance in 
the two channel Kondo model as it is screened successively by its two channels. Our model is based 
upon large-N approximation in which the localized spin degrees of freedom are described using the 
Schwinger boson formalism. 



PACS numbers: 

I. INTRODUCTION 

The Kondo effect arises from the resonant spin-flip 
scattering of electrons off localized magnetic impurities, 
which manifest themselves through the anomalous trans- 
port and thermodynamic properties of dilute magnetic 
alloys^. A modern context for the physics of local mo- 
ments is found within quantum dots 2 - in where the Kondo 
effect is manifested as a zero-bias anomaly of the dif- 
ferential conductance. The high degree of tunability in 
quantum dots allows for realization of a number of in- 
teresting variants of the Kondo effect, including both 
the "over-screened" and "under-screened" Kondo mod- 
els. For example, the Kondo effect has been observed in 
spin-1 quantum dots containing an even number of elec- 
trons, by tuning the dot to a singlet-triplet degeneracy 
point by applying a magnetic field". 

The physics of a single spin 1/2 magnetic impurity, 
with no orbital degrees of freedom is described by the 
canonical one-channel Kondo model, where the loan spin 
is fully screened by the conduction electrons at low tem- 
peratures, forming a paramagnetic ground-state with 
excitations described by Landau Fermi liquid theory. 
Nozieres and Blandin pointed out£ that in more realistic 
magnetic ions the orbital structure of the local impurity 
electron needs to be taken into account (I ^= 0), giving 
rise to a much richer class of spin-screening phenomena. 
For a spin 5 with unquenched orbital angular momen- 
tum there are ni = 21 + 1 screening channels. Depend- 
ing on whether ni is greater than, equal to, or smaller 
than the number n e = 25 of Hund's coupled electrons 
in the local moment, the spin is said to be overscreened 
(21 + 1 > 25), fully screened (21 + 1 = 25) or under- 
screened (21 + 1 < 25). In both the underscreened and 
overscreened cases, the excitations of the ground-state 
are no longer described by a Fermi liquid, and a variety of 



exotic non- Fermi liquid (NFL) phenomena can develop, 
as summarized in Table I. Unfortunately, most transi- 
tion metal ions are described by the fully screened case, 
where 25 = 21 + 1, so the underscreened and overscreened 
Kondo effect is not readily observed in bulk materials. 
For this reason, quantum dots provide an important al- 
ternative milieu for the study of unconventional Kondo 
screening phenomena. 

Most recent interest has focused on the case of the 
overscreened multichannel Kondo model. From detailed 
Bethe ansata^ and conformal field theory solutions^ the 
low-energy physics of the overscreened 5 = 1/2 Kondo ef- 
fect is known to involve a non-trivial fixed point, with sin- 

4 

gular power-law behavior of the specific heat C ~ T™!+ 2 
and magnetic susceptibility. (Logarithmic corrections to 
the specific heat appear in the two-channel case). A pos- 
sible realization of this kind of behavior in the heavy elec- 
tron compound U Be\z has been proposed by Cox, where 
quadrupole fluctuations of singlet uranium ions may as- 
sume the role of spin fluctuations in the conventional 
Kondo effect^. 

In point of fact, a spin-1/2 quantum dot coupled to two 
leads does not provide a realization of the two-channel 
Kondo effect, because the localized spin couples to a spe- 
cific combination of the lead electrons forming a one- 
channel Kondo models. A proposal to overcome this 
problem through the coupling of a quantum dot to a 
strongly interacting "quantum box" has been advanced 
by Oreg and Goldhaber-Gordonii, and the first, prelimi- 
nary indications of the the two-channel Kondo effect have 
been recently reported^ 2 -. 

By contrast, realization of the underscreened Kondo ef- 
fect requires a quantum dot with a higher spin 5 > 1/2. 
In most quantum dots, electrons fill the dot in accordance 
with the Pauli principle, forming a singlet state when the 
number of electrons is even, and for the most part, zero- 
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TABLE I: Ground state properties of the Kondo model with impurity spin S and number of channels n;. C v is the specific heat, 
and x the magnetic susceptibility. In the over-screened Kondo effect, electron-scattering remains inelastic in the ground-state. 
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bias anomalies are observed in dots with odd numbers 
of electron o 13 i 14 . A sequence of more recent experiments 
have shown that zero-bias anomalies associated with a 
Kondo effect can also occur in quantum dots with even 
occupancies, where Hund's coupling between the elec- 
trons lead to novel degeneracies, either through the for- 
mation of higher spin states, or through the accidental 
degeneracy of singlet and triplet states. The observation 
of a zero-bias anomaly in even integer quantum dots was 
first reported by Schmid et alt^ who tentatively identi- 
fied the phenomenon with a triplet ground state of the 
quantum dot. Sasaki et a& later discovered a zero-bias 
anomaly in even electron quantum dots that are tuned to 
the degeneracy point between singlet and triplet states. 
Most recently, Kogan et have shown that the singlet 
- triplet excitation energy in lateral quantum dots can 
be tuned by the gate voltage, explicitly demonstrating 
that the zero bias anomaly develops once the triplet state 
drops below the singlet configuration. 

Pustilnik and Glazmar>i£ have shown that a spin-1 
quantum dot is associated with two screening channels 
which fully screen the local moment at the lowest tem- 
peratures. The two channels arise because the two elec- 
trons in a spin-1 dot occupy two separate orbitals which 
couple differently to the leads. At low temperatures the 
conductance of the Fermi liquid that emerges is governed 
by the interference between the two channels 

G = 2^-sin 2 (<5 1 -<5 2 ) 
h 

where b\ and 82 are the scattering phase shifts of the two 
screening channels^. According to this line of reasoning, 
the development of a unitary phase shift in each channel, 
Si = 62 = 7r/2 leads to a complete suppression of the zero 
bias anomaly in a triplet quantum dofe^. 

The two separate screening channels of a spin-1 quan- 
tum dot are associated with two different antiferromag- 
netic coupling constants J\ (A = 1,2). Since the Kondo 
temperature depends exponentially on the coupling con- 
stants, Tk\ = D^/J\pexp(—j^) (p is the density of 
states and D is the band width), a rather modest dif- 
ference in coupling constants can produce Kondo tem- 
peratures separated by decades of temperature, Tki >> 



Tk2- Over the intervening temperature range, defined 
as ln(T KX ) > ln(T) > ln(T K2 ), the physics of the dot 
is expected to be dominated by the underscreened fixed 
point, or one-channel spin-1 Kondo model. In the un- 
derscreened spin-1 Kondo, the residual spin- 1/2 moment 
is ferromagnetically coupled to leads, with a coupling 
that scales asymptotically towards zerrA Recent stud- 
ies of this problem have argued that the ground-state 
which develops is a "singular Fermi liquid" , in which 
the electrons are elastically scattered with unitary phase 
shift, but in which the logarithmically decaying coupling 
to the partially screened local moment leads to a sin- 
gular energy dependence in the scattering phase shift 
and a divergence in the resulting quasiparticle density 
of state a 18 ' 19 ! 20 . This singular behavior, has a number 
of observable consequences for the conductance which if 
confirmed, will provide a first experimental realization of 
singular Fermi liquid behavior of an underscreened local 
moment. 

In this paper we present an analysis of the transport 
properties of a spin-1 quantum dot in two cases: first 
for one screening channel and then for two active screen- 
ing channels. The initial part of our paper provides a 
detailed treatment of the singular conductance proper- 
ties previously predicted by us in a short paper on this 
topic 21 . The second part of the paper extends this ear- 
lier work to take account of the second-channel screening 
processes and the interference it gives rise to. One of the 
main results of this new analysis, is the development for 
an approximate expression for the temperature depen- 
dence of the linear conductance. We find that the linear 
conductance can be divided into a coherent part, and an 
incoherent part, 

G = G co h + Gi nc 

where 

G coh = ^(2aP) 2 J oIlo (-^) Wi(w)-r a («))f 

defines the coherent conductance through the two- 
channel and 
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x ^[r^M + ^V)] (1) 

the additional contribution due to incoherent transport, 
where "Ji" and "T2" refer to the scattering t-matrices in 
the two screening channels of the quantum dot, a and f3 
define the relative amplitudes of channel f in the left and 
right channels respectively and 

r A - x =ImT A H-^|T A H| 2 . (A =1,2) (2) 

are the inelastic scattering rates in each channel, defined 
by the deviation of the imaginary part of the t-matrix 
from the value expected from the optical theorem. 

II. CONSTRUCTION OF THE HAMILTONIAN 

Our initial model Hamiltonian for a spin-1 quantum 
dot coupled to two leads can be divided up into a Hamil- 
tonian for the leads {Hl), the quantum dot (Hd) and the 
coupling between the leads and the dot (He) as follows 

H = H L +Hc + H D (3) 

where 

H L = ^c\ ia c kla (4) 

fctr 

y=L,R 

describes the electrons in the left and right leads, 

H c = h ]T(^L^ + Rc -) +*2 E^L d 2 CT + H.c), (5) 

k k 

describes the hybridization in two channels between the 
lead electrons and the quantum dot through the contacts. 
Here d\ a and d\ a create electrons in two orthogonal states 
of the quantum dot and 

V>L = aC kLo+P C kRo 

tp k(T = -l3c\ LlJ + acl Ra , (6) 

create the linear combinations of conduction electrons 
that hybridize with these two orthogonal scattering chan- 
nels. The Hamiltonian of the spin-1 quantum dot can be 
written 

H D = -Exndi - E 2 n d2 + Uin d i^n dll + U 2 n d2 ^n d2l 
+ Ui 2 n dl n d2 - J [Si ■ S 2 }. (7) 

Here — E\ (A = 1,2) are the energies of the two one- 
particle states and n d \ = J2\ d\ a d\ a (A = 1,2) are the 
occupancy of the two channels. U\, U 2 and Ui 2 are the 
intra and inter-channel Coulomb interactions and lastly 
Sa = \d\ a S a f}d\f} are the spin operators for the two 
channels and J is the (ferromagnetic) Hund's exchange 
coupling between the spins. 

In general, we are interested in the case where the one- 
particle energies — E± and — E 2 are negative, but Ux, U 2 
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FIG. 1: Energy diagram of a spin-1 quantum dot with strong 
Hund's coupling. 

are large quantities that restrict double occupancy of ei- 
ther d-state. Finally, the Hund's interaction J is assumed 
to be large restricting the doubly occupied state to be a 
spin-1 triplet state. In this situation, the quantum dot is 
restricted to the states 

\d°) = |0), 
Id 1 :1a) = dl\d°), \d l :2a)=dl\d°), 

r 4 T 4 T [d°>, (m = i) 

\d 2 :M) = I ^(4|4 T +4 T 4l)K), (^ = 0)(8) 

l4|4|M°): (M = -l) 

The energy of the d 2 triplet states is E\ + E 2 — J/4 + U\ 2 
(see Fig.l). 

We can project the Hamiltonian onto this Hilbert space 
by combining the d-electron operators with a generalized 
Gutzwiller projection operator as follows, 

4« - Prdl = Xl (9) 

where 

Pt = PgiP G 2 (l + Si • S 2 - (10) 

and 

P GX = (l-n m n dXl ), (A- 1,2) (11) 

is the Gutzwiller projector for channel one and two. Here, 
we have used Hubbard's "X" notation to describe the 
projected d-electron fields. With this notation, we can 
write our projected Hamiltonian as 

ka ka 

+ h J2^ta X f + h - C -) + *2 J2(vtX 2u + h.C.) 

k k 

- El ^ X^Xig - E2 2J X 2a X 2<y 

+ (Uu - n dl n d2 (12) 

The main interactions in this Hamiltonian are hidden in 
the constraints. In the Kondo limit, where this Hamilto- 
nian is dominated by excitations between the d 2 and d 1 
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state, the final interaction term has the effect of replacing 
the excitation energy 

-E x -» E(d 2 ) - E{d l ) = -E x -j + Ux2 

With this observation, we shall eliminate the final inter- 
action from the Hamiltonian, absorbing its effect into 
the above redefinition of E\. (Another way is to set 
012 = J/4). 

In the extreme Kondo limit where the dot is predomi- 
nantly in the \d 2 ) state, the d 2 <-> d 1 charge fluctuations 
can be eliminated from this Hamiltonian via a Schrieffer 
Wolff transformation to give 

H = H L + J X ^ i>i a Vaf3lpk> 13 ■ S 

+J 2 E <t>\ a a a f)<t>k>i3 ■ S (13) 

k,k'a(3 



where 



J 2 = 



(t2? 



i + E 2 -U 12 ' 



(14) 



are the two antiferromagnetic coupling constants, corre- 
sponding to two orthogonal channels. 



= -j=b}bl\d u } = -j=\d':m = 0).(W) 
which consistently normalizes the matrix element 



< d 2 : m = 0|Zj T |d x : 1 |) = -^=. (19) 
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The physical Hilbert space is defined by those states 
where 



nb + n x i + n X 2 = 2S 



(20) 



A. Large-N expansion 

To develop a controlled many body treatment of JT2|) 
we shall employ a large- N expansion, extending the num- 
ber of spin components a from two to N. For the work 
in this paper, we are interested in the approach to the 
Kondo limit of this problem, where the amplitude of the 
valence fluctuations are small, and in this limit, we shall 
replace y/rib — * V / 25 in (fTS)) the above expression. With 
this device, 



^=btxx, (A = 1,2) 



III. SCHWINGER BOSON REPRESENTATION 

A very convenient way to treat the Gutzwiller opera- 
tors in this Hamiltonian is to use a combination of slave 
fcrmions and Schwinger bosons. If we describe the empty 

state \d°) = X2X1IO) as a P arr 01 two slave fermions, then 
the singly occupied states are written 

\d l :la)=bixi\d°) = -btxt\0), 

Id 1 : 2<r> = &t X2 |d°) = 6t x t| 0)) (15) 

while the doubly occupied states are given by 

\d 2 : M) = {^(^t) 2 |0), b\b\\0), ^(^) 2 |0),} 
In this representation, the Hubbard operators are written 
Xl = -±=bixu Xl = -^btx2, (16) 



The prefactors ^= (n& is the number of Schwinger 

bosons) are needed to guarantee the normalization of the 
wavefunction. Thus 

\d 2 : M = 1) = X^X^O) = ^(b\) 2 \0) (17) 
which is the right normalization. Similarly 

p T d\/ 2l \d°) = x^xlid ) 



where N is the spin degeneracy and t\ = t\/ \/2S/N . 
The model Hamiltonian is then 



ka 



-h.c.) + ^=J2&txW + h.c.) 



Eixlxi + E 2 xlx2- 



(21) 



where it is understood that as N becomes large, ii,2 is 
kept fixed. There is one additional trick required to pre- 
serve a finite scattering phase shift as N — > oo. For this 
task, we introduce K — kN bosonic "replicas" , where k is 
maintained fixed as N — > oo. With this device, in strong 
coupling K bosons bind into the Kondo singlet and we 
can produce a large N mean field theory in which the 
scattering phase shift is S = irk, with the qualitatively 
correct logarithmic energy dependences 20 . The Hamilto- 
nian used in the large N expansion is then 



ka 



ka 



ka/x 
ka/x 



E 
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+ X(n b + n x i + n x2 - 2SK) 



(22) density of states, N(u>) 



where the sum over /i runs from 1 to kN, while the last 
term imposes the generalized constraint 



rib + n x i + n X 2 = 2SK. 



(23) 



Notice how the large N model Hamiltonian (|22|) is ob- 
tained by replacing 



N 



X 



b Ux^- 



(24) 



in (TT2"]) . Section IIVI examines in the underscreened one- 
channel regime, where t 2 — 0. Section [V] proceeds with 
a discussion of the full Hamiltonian. 



IV. ONE-CHANNEL SPIN-1 QUANTUM DOT: 
UNDERSCREENED KONDO EFFECT 

A. The singular Fermi liquid 

In this section we discuss the limit where the Kondo 
temperature of the second-channel Tki is negligibly 
small, so that the physics of the quantum dot is governed 
by an underscreened Kondo effect. We present a large N 
approximation to the physics of this underscreened limit. 

The Hamiltonian that governs this behavior is derived 
from the Hamiltonian (|22)l with t 2 — and E 2 = 0. Much 
is known about the equilibrium physics of this model. At 
low temperatures, the spin is partially screened from spin 
S to spin S — (1/2). The residual moment is ferromag- 
netically coupled to the conduction sea, with a residual 
coupling that slowly flows to weak coupling according to 



tK, which means 



1 d8(uj) 

7T dui \uj\ ' 

that we can not associate this state with a bona-fide Lan- 
dau Fermi liquid. For this reason, the ground-state of 
the underscreened Kondo model has recently been called 
a "singular Fermi liquid"—. 

These singular features of the underscreened Kondo 
effect are expected to manifest themselves in the prop- 
erties of a triplet quantum dot, in the range where 
Tk2 << T << Tki- In this regime, the physics is domi- 
nated by the underscreened fixed point. For example, as 
a function of magnetic field, we expect the conductance 
to follow the simple relation 



G{B) 



dl 
~dV 



• sin 2 6(B) 



so at low fields Tki « HbB « Tki, the differential 
conductance will have the form 



1 - 



1 
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Notice that the field derivative of the conductance di- 
verges as 1/B 



dG 
dB 



If! 
B~h 



8lr 



Tki 



at low fields. The prediction of the finite temperature, 
and finite voltage conductance can not be made exactly, 
however we expect the above form to hold, for the differ- 
ential conductance at finite temperature or voltage, with 
an appropriate replacement of cut-offs, namely 




Jp(A) 



where A ~ max(T, iibB) is the characteristic cut-off en- 
ergy scale, provided in equilibrium, by the temperature 
or magnetic field. At low energies and temperatures, the 
partially screened magnetic moment scatters electrons 
elastically, with a unitary phase shift, however the cou- 
pling to the residual spin (S — ^) gives rise to a singular 
energy dependence of the scattering phase shift. The low 
energy scattering phase shift can be deduced from scaling 
theory to have the asymptotic form 



i ± 



(S- 



]n(T K1 /u>) 



for the "up" and "down" spin channels. The logarithmic 
term on the right hand side is produced by the residual 
coupling between the electrons and the partially screened 
moment. While the electrons at the Fermi energy scatter 
elastically off the local moment with unitary scattering 
phase shift, as in a Fermi liquid, the logarithmically sin- 
gular dependence of the phase shift leads to a divergent 



G(V,T)~- 



1 - 



f 



16 



( T K1 \ 
^max(T,eV) J 



so that over the temperature range dominated by the un- 
derscreened Kondo fixed point, the conductance will be a 
monotonically increasing function of decreasing temper- 
ature and voltage. Over this range, dG/dT ~ 1/T and 
dG/dV ~ 1/V will be divergent functions of temperature 
and voltage respectively. 



B. Dyson equations and one-channel t-matrix 

We now analyze the low energy behavior of the t- 
matrix in the large-N approach. Our development be- 
gins with a derivation of the Dyson equations for the 
Green's function of the conduction 'i/j/cr-fermions (G) that 
are coupled to the dot and the Green's function of aux- 
iliary x— holons (J) that describe the partial screening. 
In this large-N approach, we will show that the Green's 
function of the Schwinger bosons B remains unrenormal- 
ized. 
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In general, we will be dealing with a voltage-biased 
quantum dot. The Hamiltonian for the leads is given by 

H a = J2(tk - fJ, 7 )cl ia c kjcr , {j = R,L) 

where /xj, — and pn = —^f. The Green's function 
of the conduction electrons is given by 

9j(t - t 1 ) = -i^(T CK c^(t)4 7a (t')) (25) 
k 

where ck is the Keldysh contour. Following the stan- 
dard procedure, we write this Green's function using the 
Larkin-Ovchinnikov representation, in terms of the ad- 
vanced (g A ), retarded (g R ) and Keldysh Green's func- 
tions g K a,s follows 



and 



5 7 









9* 









9$_ 





—inp —i2irp(l — 2/ 7 ) 
iirp 



where p is the electron density of states and 

f 7 (e) = f(e-^), (j = R,L) 

describes the Fermi function in the left and right-hand 
leads. The complete Green's function describing the two 
leads is then a four-dimensional, block-diagonal matrix 



9le 




(26) 



To obtain the corresponding Green's function in the 
channel basis, we carry out a unitary transformation, 
writing 



9ch = 



gipip 9^<p 
9<pi(> 9<pp 



— RgieadR 1 



define weighted-averages of the distribution functions in 
the two leads. The subscripts (0) are included, to delin- 
eate these functions from the renormalized local conduc- 
tion electron distribution functions that can, in principle 
develop when the leads are coupled. 

There are no retarded or advanced components to the 
off-diagonal block matrices g c h, but the inter-channel 
Keldysh Green's function does become finite when there 
is a voltage between the two leads, 



9U 




where 



4» = = (to*p)<*0(fR(u>) - /lH), (30) 

and in this way, a finite voltage mixes the two channels. 

We shall return to these general expressions later. 
However, for the one-channel problem, only the if; elec- 
trons couple to the quantum dot, and so our interest now 
focuses on g^. 

The Dyson equations for the propagators of conduction 
^-electrons and slave x-fermions (Fig. 2) are 



G- 1 
J- 1 



9^\ — ^1 



where 



-iirp 





n, 



1 



(31) 
(32) 



(33) 



is the bare inverse conduction propagator and 



where 



R 




(27) 



transforms between the two bases. The block-diagonal 
components of the matrix g c h are now 



gi>i> 



and 



where 



fi \u>)=a?f L (u;)+lPf R (w) 



-iirp -2mp{\ - 2/< 0) ) 
iirp 



-iirp -2iTrp{l-2f!p ) ) 
iirp 



(28) 



(29) 



Jo 1 



(u-E d -\ + i8) 

(w - E d - A - i8) 



(34) 



is the inverse propagator of the \ fermion. Notice how the 
Keldysh component of = ^4^(1 -2 f ^ (u)) is ^itc, 

whereas the Keldysh component of J^ 1 is infinitesimal, 
and has been dropped from the above equations. 

The self-energies are generated in our large-N mean- 
field theory within the non-crossing approximation 
(NCA). The diagrammatic Dyson equations are shown 
on Fig.l. The boson is our large-N approach is approxi- 
mated as a sharp excitation, with an average occupancy 
(n bo> ) = n B = (2SK)/(NK) (see Appendix I). In all our 
calculations, we make use of the bare Green's function for 
the Schwinger bosons 



B, = 



B£ 



(35) 



7 



ca-e 



The self-energies for the conduction electrons are 



0) 

ca-e 



+ 



0) 



FIG. 2: The non-crossing approximation for the self-energies 
of conduction electrons and \ fermions. The solid line denotes 
the Larkin Ovchinnikov matrix propagator for the conduction 
electrons. The dashed line denotes the corresponding Green's 
function of the auxiliary (\) fermions and the wavy line is the 
bosonic propagator. Thin lines denote the bare propagator 
and full lines the dressed propagator. Each vertex corresponds 
to the factor i-^=- 



where 



B, 



R,A 



v — X±i5 



(36) 



where h B = 1 + 2ns (A). 

/From the Dyson equations we obtain sets of self- 
consistent integral equations for the retarded, advanced 
and Keldysh self-energies (Appendix I) . The self-energies 
for the slave (\) fermion are 



Il R (e) = -t 2 n B (X)G A (X-e) 



J-oc T 



-JmG R (uj) 



uj — A + e + id 
U K (e) = -2it 2 ImG R (X- e)[n B (X) - U(X- e) 

-- 2n B (\)U(\~e)}. (37) 

The advanced self-energy is determined from the complex 
conjugate of the retarded self-energy H A (ui) — U r (uj)*. 
In obtaining these results, we have assumed that the con- 
duction electron Keldysh Green's function is determined 
by the relation 

G k (uj) = (G r (uj) - G A (uj))h 4 ,(uj) 

where a priori, h^(u>) = 1 — 2f^(u>) no longer equal to its 
equilibrium distribution. 

The ratio of the Keldysh to the retarded self-energies 
self-consistently determines the x fermion distribution 
functions 



n 



A' 



h^h B — 1 



n R -if 4 h B - 

Rearranging this expression gives 

h x h B — 1 



h^(uj) 



h B — h x 



(38) 



(39) 



£ R (e) = - t ^n B (X)J A (X-c) 



k r°° dio 



,K 



uj — A + e + id 



-ImJ R {uo) 



£ A » = -2it 2 -ImJ R (X-e)[n B (X)-f x (X-e) 
- 2n B (X)UX-e)\. 



(40) 



and again, E A (w) = T, r (uj)*. The ratio of Keldysh to 
the retarded/advanced Green's function is given by 



/ty (uj) 



-*K 



Z R - ^ A h B - h x 



(41) 



but this recovers exactly the result obtained in (|39[) . 
showing that detailed balance is satisfied. However, we 
can't choose any distribution function hf. If we go back 
to the original Dyson equation for the conduction elec- 
tron Green's function (|3"Tj) , the Keldysh component of 
the Green's function is given by 

G K = G R [X K - (g^) K ]G A 

But the distribution function associated with the self- 
energy is determined by "E K = h^(Jl R — T, A ) whereas 
the distribution function associated with [g^] K = 

2/(—inp)(l — 2/1 ) is the bare distribution function 



,(0) 



2/r 



In this way, we see that our original 

(0) 



assumption G K = (G R — G A )h^ requires that — K. 
In other words, in large N limit of the single-channel 
quantum dot electron distribution function is unaffected 
by the coupling to the dot. 

We can summarize the results of our calculation of the 
Keldysh self-energies by providing the distribution func- 
tions that they generate. The distribution function of 
the auxiliary fermion is simply determined by the rela- 
tionship 



n b [l - fa(u))] 
n b + J^(uj) 



(42) 



where ra& = l/(e' 3A + 1) determines A. This relationship 
can be simply understood as the result of detailed balance 
between rate of the decay processes c — > b + \ and b + 
X —> c, and it reverts to the equilibrium Fermi Dirac 
distribution in the limit V — > 0. 

We have solved the equations (|3"7) and (|4T]|) numerically 
to obtain the energy-dependent i-matrix, given by 



i^(u>) = S(w)(l - g^(oj)Y,(u>)y 



(43) 



Actually, the Keldysh part of this quantity is determined 
simply from the relation = (t R — t A )h^\ and the ad- 
vanced part of the electron t-matrix can be simply writ- 
ten 



(1 - iirpZ A (uj)) 



(44) 
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eV/T K 

FIG. 3: Imaginary part of the t-matrix for a variety of voltages 

for K/N = 0.4. As the voltage is increased, the singular FIG. 5: Voltage-dependent conductance G(V) = I(V)/V for 
central peak splits into two components. the case K/N — 0.4. Insert: derivative of G(V)showing 1/V 

divergence. 



g(T) 

Hu 




Tk/T 



T/T K 



FIG. 4: Temperature dependence of the differential conduc- 
tance, normalized with respect to gu — N ^-sin 2 (irK/N) for 
the representative case K/N — 0.4. Insert shows the 1/T 
divergence of the derivative dg/dT. 



The imaginary part of t A {u>) determines the electronic 
density of states of the resonantly screened d-electrons. 
The voltage dependence of this quantity (at T = 0) is 
shown in Fig. 3. At zero voltage, the t-matrix contains a 
logarithmic divergence, which splits into two peaks at a 
finite voltage. In this large N limit, the split Kondo res- 
onance retains its singular structure. This is an artifact 
of taking a large N limit in which the Schwinger bosons 
are preserve their sharp spectral structure. 



C. Conductance 

The conductance G = 31/ dV is defined by the tem- 
perature and voltage dependent current 2 ^ 



I(V,T) 



2h 



Np 



duj 



[/iH-ZflMjlm^H (45) 



where the advanced t-matrix is given by (I44|) . 

In Fig. 4, we show the computed temperature de- 
pendent conductance. The temperature-dependent 
deviations from unitary conductance are determined 



by the logarithmic singularity in the phase shift, 
and in our calculation, these are proportional to 
1/ ln(TV/max(e V, T)). There is a technical point here 
that needs some discussion. In the Schwinger boson ap- 
proach, the number of bound bosons in the Kondo singlet 
does not exceed N/2, which has the largest bound-state 
energy, so the region K > N/2, does not propertly de- 
scribe an underscreened Kondo model. Consequently, we 
are limited to static phase shifts S = n(K/N) < tt/2, so 
the particle-hole symmetric case 5 = tt/2 lies beyond the 
reach of this approach. This will leads to some important 
differences between the results of the current calculation 
and those expected in experiment. In general our solu- 
tion does indeed capture the singular low energy behavior 
of the phase shift 



irk 



cons 



(46) 



Now the conductance G depends on sin 2 S(V, T), so that 
for k ^ 1/2 



G ~ sin 2 5 = sin 2 (7rfc) 



cons 
In [2H 



(47) 



which is the form that our solutions follow. However, in 
the special case k — 1/2, the coefficient of leading log in 
conductance vanishes, so the conductance involves square 
of the logarithm 

G - sin 2 (7r/2) 



In 2 [^] 

In both cases, the temperature derivative of the conduc- 
tance has the singular divergence dG/dT ~ 1/T, a fea- 
ture which can be tested experimentally. 

Finally, Fig. 5 shows the voltage dependence of the con- 
ductance, which has a similar logarithmic singularity at 
low voltage. The insert shows d 2 I/dV 2 versus 1/eV. 
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V. TWO-CHANNEL SPIN-1 QUANTUM DOT: this section we consider the system, described by the 
TWO-STAGE KONDO EFFECT Hamiltonian d5TJ and finite voltage. We show that at 

finite voltage the current has three contributions: contri- 
We start with a derivation of the general current ex- bution from each channel and the interference term, 
pression for the case of a two-channel quantum dot. In 

I 



A. Derivation of the two-channel current expression 

The current from the dot into the left, and right hand leads is defined from 

Il,r = -e(N L , R ) = --([H,N LtR \). 



(48) 



The non-zero contribution to the current derives from the commutator with the hybridization Hq in equation ([5]) . In 
equilibrium I R = —II = I- The current into the right-hand lead is given by 



In 



[PH X L c kRv) + at 2 {Xl a c kRa ) - H.c 

ka 

-J E* 1 [a/3((*iW> - H.c.) +[3 2 ((xljj ka ) - H.c.) 

ka 

f£*3 [ a /3«xU feCT ) - H.c.) + a 2 ((xlvk*) - H.c. 



(49) 



Now the expectation of the occupancy of the dot in each channel is a constant in the steady state, so that 

-{n dX )= % -{[H c ,n A> ]) = Q, (A = 1,2) 



Carrying out the commutator, we obtain 







-H.c 



E K 



(fika) - H.C 



ka ka 

From this consideration, we see that the second terms in equation ([49]) identically vanish, so that the current 

I = Ir = - apJ2 + H.c.) + t 2 ((xlij ka ) + H.i 

ka 



(50) 



(51) 



A similar procedure for Jj, confirms that Jl = —I. Note, that in contrast to the Hamiltonian, the current is determined 
by the overlaps between the ^-field with the cfe-electron an d the (g- field with the d\ -electron on the dot. 



Combining the Dyson equations for the Green's func- 
rectly related to the equal-time Keldysh Green functions, tions G± v and with the expressions for bare Green's 

functions from Section IVB (for details see Appendix II) 
we obtain the general current expression in case of two 
channels 



The expectation values that enter into current are di- 

ions, 

(52) 



ka 



and 



ka 



enabling us to write the current in the form 
+t 2 (G^ 2 H-G^M)] 



T Ne R 



(53) 



(54) 



I = N{I V + I^ + I int ), 



(55) 



where 



2e, 



{aPYtiP / Mfdto) - f R (u))ImD& 



h = \{^fi\p J du(J L (u) - f R ((4))ImD&, (56) 

where D i3 {t - t') = (T cK X ia (t)xj a (t')) is the Green's 
function of the electrons on the dot. 
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The interference between the channels is induced by 
the terms 



hut = —^.aptxhip J du{D& + D& 

+ (1 - 2/3 2 /lH - 2a 2 f R (u)){D^ - D*) 
+ (1 - 2a 2 /lM - 2(3 2 f R (u;))(Df 2 - £>£)} 



(a) 



(57) 




Notice that in the one-channel limit, (f 2 = 0), the total 
current reduces to 



I = NI i> = ^-[ duj[f L (uj) - f R (u)]~ImD&(L 



(58) 



recovering the Meir-Wingreen^ 2 - expression for the cur- 
rent in a one channel quantum dot, with T = 
r L lV(r L + T R ), where r L = 2npaHf and T R = 
2np(3 2 t{. 



B. Approximate expression for the conductance 

In order to calculate the conductance, we need the 
Dyson equations for the Green's functions, defined by 
Hamiltonian (|2"2"|) . The Green's functions are defined in 
the same way as for the one-channel case and the Dyson 
equations (Fig 6) are 



Gipij} G 






] _ 


| E^ Tj^ap \ 










Here 



(io _1 )n o 

Uo'h 



(J^' A = LJ-E i -X±iS 



0. 



(59) 

rin n 12 \ 
n 2 i n 22 J 

(60) 



(61) 



Unlike the single channel case, the solution of these equa- 
tions at finite voltage is complicated by the need to self- 
consistently compute the distribution functions. This is 
because the hybridization no longer commutes with the 
Keldysh Green's functions of the conduction sea. In the 
following we limit our calculation to the linear response 
at zero voltage, which can be written in terms of the 
equilibrium V — Green's functions. 

In equilibrium at V = 0, the Green's functions are 
channel-diagonal 





~ 9ipip 






A-l 

<P<P 


9 '<p<p 








= U°n)- 


1 flu, 




J 22 


= 142)- 


1 -n 22 . 


(62) 




(b) 



2 2 



1 1 



FIG. 6: Dyson equations for the conduction electron propa- 
gators (a) and the slave fermion propagators (b) for the two- 
channel quantum dot. 



In this case the single-channel contributions to the cur- 
rent (|56p can be expressed via "channel" t-matrices 

h = ^(«/3) 2 / MfUoj) - f R {u))ImT^{w) (63) 
I v = ^(uP? I Mh(u) ~ f R (u))ImT v (u), (64) 



where = E^(l - iTrpE^)" 1 , and T v = E^(l - 
inpT,^)" 1 are the advanced t-matrices. In linear re- 
sponse the single-channel t-matrices are proportional to 
the diagonal Green's function of the electrons on the dot 
holds a simple relation between channel Green's func- 
tions and the Green's functions of the physical electrons 
on the dot 



T^u) = (ti) a Dft(w), 



(65) 



where 



D$ x ,(uj)=i ({X Xa (t), Xl(0)})e^dt (66) 



and X L = vIsS^i^^Xm- In the lar S e N limit ' 
the t-matrices for each channel completely decouple in 
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equilibrium, and are given by the scaling form 



T x (u,T) = T( 



T 



Tkx T, 



KX 



) 



(67) 



where Tkx = 1, 2) are the Kondo temperatures in the 
two channels. 

The interference contributions to the current pose a 
greater challenge and we have been unable to treat these 
terms without making some additional degree of approx- 
imation. Without loss of generality, the Keldysh compo- 
nents of the dot Green's functions can be written 



d k = D it(_ D 



-UK D A 



We can further simplify the full expression for the lin- 
ear response current by taking advantage of the optical 
theorem. In the case where the scattering off the im- 
purity is purely elastic, the channel diagonal t-matrices 
satisfy an optical theorem ImT\(ui) = irp\T\(ui)\ 2 . The 
deviation from the optical theorem describes a decoher- 
ence scattering rate-' as follows 



ImT A H-^|T A (w)| 5 



(73) 



In this form, the combined interference and channel di- 
agonal currents can be written in the simple form 



G = G C 



Gi- 



la linear response, the effect of the voltage on will appear 
through the voltage dependence of (D^ 1 ) K , so we can 
replace the retarded and advanced combinations in the 
above by their channel diagonal components, 



where 

Ne 2 

G coh = — (2a(3) 2 I dio 



dm 

du) 



Df x {-D 



-UK tjA 



(68) 



In general, the Keldysh components of and D^ 1 
contain a hybridization component and a term derived 
from interactions 



defines the coherent conductance through the two- 
channel and 



Gin,c. — 



(2a/?) 2 



duj — 



du> 



X TTpfTj (w) + r 2 



(74) 



{-D^) K = t l9 ^h - (f x t a )E 



(69) 



The second terms are a kind of vertex correction to the 
current operator (see Fig. 10 (a)). These terms describe 
the inelastic corrections to the interference contribution 
to the current. Our approximation entails neglecting 
these vertex corrections, so that (see Appendix III) 



= titzDfig^Dyx 



D K 



ht 2 D*g 



K jjA 



(70) 



When the same vertex corrections are ignored, the re- 
tarded and advanced Green's functions are channel- 

-,R,A 



diagonal, D^ A 



= o. 



With these approximations, the contribution to the 
current due to the interference between the two channels 
becomes 



lint 



(71) 



or, 



2e. 



i( M 2 V| du(f L (uj) - Jr{uj)) 

[T^H^H+Tx^^H]. (72) 

This expression brings out the interference character of 
the term, and combined with (|63| describes the linear 
response current flow, purely in terms of the channel di- 
agonal, equilibrium t-matrics. 



the additional contribution due to incoherent transport. 
In the limits where the scattering t-matrix is purely elas- 
tic, 



Tlf) 



so that the inelastic scattering rate r A 1 = so at low 
temperatures 



Gr 



Ne 2 
Ne 2 



(2a/3) 2 M 2 |T 1 (0)-T 2 (0)| ; 
(2a/3) 2 sin a (*i-*0 



(75) 



recovering the result expected from Landauer theory. 

In the large N limit it is straightforward to compute 
this approximate expression for the conductance. The 
scaling form (|67|) can be used to rescale the t-matrix from 
one channel to the other. Thus if Tkx/Tk2 = ol, then 
from (f6"T)) . 

T 2 (w,T) = Ti(au,aT). 

The behavior of the linear conductance in case of the 
two-channel quantum dot for different ratios of the two 
Kondo temperatures is shown on Fig. 7. When the two 
Kondo temperatures are equal to each other, the conduc- 
tance is completely suppressed due to destructive inter- 
channel interference. When Tkx and Tki are widely sep- 
arated, at low temperatures conductance is suppressed as 
expected, and at high temperatures it develops a char- 
acteristic hump-structure due to the Kondo resonance in 
one of the channels. 
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-3-2-10 1 2 

log 10 (T/T K ) 



FIG. 7: Linear conductance of the two-channel quantum dot 
normalized with respect to gu — N ^- sin 2 (ttK /N) , calculated 
for K/N = 0.4 and a = j3 = 1/V2, using equations (f74|) and 
(|74|) . for three ratios of the two Kondo temperatures: a — 
Tki/Tk2 = 10 (short-dashed curve), a = Tki/T K 2 = 100 
(long-dashed curve) and a = Tki/Tk2 = 1000 (solid curve). 

VI. CONCLUSION AND DISCUSSION. 

This article has considered the physics of a quantum 
dot in which Hund's coupling between the electrons lo- 
calized within the dot gives rise to a spin-1 configuration 
on the quantum dot. Theoretically, such a quantum dot 
is expected to map onto a two-channel, spin-1 Kondo 
model in which a ir/2 phase shift develops in both scat- 
tering channels. According to a Landauer analysis of the 
resulting Fermi liquid, the transmission through the dot 
should vanish at low temperatures, giving rise to a con- 
ductance 

e 2 

G = 2— sin 2 (<Ji -S 2 ) (76) 
h 

that vanishes at low temperature o 16 ' 17 . Paradoxically, 
experiment suggests that spin-1 quantum dots do indeed 
develop a zero-bias anomal y 3 ' 4 ' 15 . It is this issue that has 
motivated the current body of work. 

Motivated by these results, the current authors 
proposed 21 - that triplet quantum dots with a large zero 
bias correspond to anisotropic spin-1 Kondo models, in 
which the Kondo temperatures of the two channels ex- 
hibit a large ratio. In this paper, we have explored the 
general problem of high-spin quantum dots, screened by 
two conduction channels, attempting to develop a frame- 
work to model the detailed conductance of such systems. 
Our proposed model involves an infinite U description 
of the electrons in the quantum dot, with an infinitely 
strong Hund's interaction that ensures the formation of 
a high-spin electron configuration in the dot. In this 
case the ground state of the dot is a spin-triplet and 
the charge fluctuations to the excited singly occupied 
states can be eliminated via Schrieffer- Wolff transforma- 
tion. The Hamiltonian which describes the spin-1 quan- 
tum dot at low temperatures is then that of a two-channel 
Kondo Hamiltonian with two different Kondo coupling 



constants, and two different Kondo temperatures for the 
two electron states inside the triplet. 

The emergence of two scattering channels in a spin-1 
quantum dot has interesting consequences: in particu- 
lar, the conductance develops a a non-monotonic depen- 
dence on magnetic field, temperature and voltage 1 ^. At 
zero temperature, destructive interference between the 
two elastic scattering channels leads to a complete sup- 
pression of conductance and the resulting ground state is 
described within a Landau Fermi liquid picture. Con- 
ductance starts gradually increasing with temperature 
and between Tki and Tki reaches it's maximum (the 
height of the maximum which is less or equal to the uni- 
tary conductance, depends on the difference between the 
two Kondo temperatures). At still higher temperatures 
T >> Tki, the conductance is ultimately suppressed to 
zero. Note that when the two Kondo temperatures are 
equal to each other, conductance remains suppressed for 
all temperatures 1 ^. 

Our analytic analysis of this non-monotonic depen- 
dence of conductance on temperature employs a a 
Schwinger boson formalism to describe the dot spin de- 
grees of freedom, using a large-N expansion to pro- 
vide an approximate treatment of the resulting many 
body physics. Although the bosons in our approach re- 
main unrenormalizcd, which means we neglect the spin- 
decoherence effects we find that the that the initial in- 
crease, and ultimate suppression of conductance are re- 
produced in our analysis. 

As part of this work, we considered the interesting 
case when one of the scattering channels is completely 
suppressed, which provides a realization of the under- 
screened Kondo effect in the spin-1 quantum dot. In this 
case at zero temperature conductance reaches its uni- 
tary limit, but in a specific singular way. Although the 
scattering shift at T = is equal to 6 = ir/2, the energy 
dependence of S contains logarithmic contributions. This 
situation is described in terms of singular Fermi liquid 19 . 
It turns out that our Schwinger boson approximation is 
suited for description of a singular behavior of conduc- 
tance at low temperatures and low voltages. 

There are several open questions, of both a theoreti- 
cal, and experimental nature that arise from this work. 
The current through the one-channel quantum dot is 
compactly related to the density of states of the dot- 
electrons via the Meir and Wingreen relation. In the 
two-channel case, these simplifications continue to oper- 
ate in the contributions to the current that are channel 
diagonal, but the off-diagonal components to the current 
require a knowledge of the Keldysh Green's functions 
and the voltage-induced components to interchannel t- 
matrix. Is the physical reason for this added complexity? 
One possibility, that we have not been able to eliminate, 
is that a more general matrix ansatz for the Keldysh 
Green's functions, of the form, such as 

G K (u>) = G R {u)F{u) - F(lo)G a {uj) 

can be used to simplify our results. 
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On the experimental front, we are still lacking hard evi- 
dence that the zero-bias anomaly in triplet dots is driven 
by an extreme ratio of Kondo temperatures. It would 
be particularly interesting if lateral quantum dots of the 
kind used in previous experiments could be adapted to 
support a variable ratio between the Kondo tempera- 
tures in the two channels. This might be done, for ex- 
ample, through the introduction of additional gates that 
change the symmetry of the quantum dot. This might 
provide the means to observe the cross-over from the sin- 
gular Fermi liquid associated with fully developed under- 
screened behavior, to the interference dominated, non- 
monotonic conductance of the fully screened quantum 
dot. 

The authors wish to thank H. Kroha, G. Schon and 
G. Zarand for discussions and comments related to this 
work. We are particularly indebted to M. Eschrig for 
many useful discussions. This research was partly sup- 
ported by DFG CFN (AP, BB) and DOE grant DE- 
FG02-00ER45790 (PC). 

Appendix I: Self-energies for the one channel Kondo 

effect. 



a) [lO 




~ 0(K/N)~ 0(1) 



b) 

./ V ~ 0(N/N)~ 0(1) 

a 




0(1/N) 



FIG. 8: NCA contributions to the self energies of a) lead 
electrons, b) auxiliary x-holes, and c) Schwinger bosons. The 
solid line denotes the Larkin Ovchinnikov matrix propagator 
for the conduction electrons. The dashed line denotes the 
corresponding Green's function of the auxiliary (x) fermions 
and the wavy line is the bosonic propagator. Each vertex 
corresponds to the factor i-?=. 



The self-energy contributions to the equations (|31j) and 
(|32|) are shown in Fig. 8. Each vertex is associated with 
a factor -^=. Internal summations over the spin index, 
or replica indices inside the loop provide a factor of iV or 
K respectively, so that the conduction and slave fermion 
self energies are of order 0(N/N) = 0(1) and 0{K/N) 
respectively, and they remain finite in the large N limit. 
By contrast, the bosonic self-energy contains no internal 
summation over internal quantum numbers, so that this 
quantity is of order 0(1/ N) and vanishes in the large N 
limit. 

The conduction electron self energy is directly related 
to the Green's function of the d-electron, written in com- 
posite form as 

-iP 



(i-0 



N 



(77) 

where T CK refers to time-ordering along the Keldysh con- 
tour and the notation A± = A\ ± A-i denotes the sym- 
metric and antisymmetric combination of fields from the 
forwards time (1) and backwards time (2) parts of the 
contour. (Sometimes A + = A c i and A_ = A q are referred 
to as the classical and quantum fields, respectively.) 

The retarded, advanced and Keldysh combinations of 
the self energy are related to 



= Ej 



E A = E_ 



(78) 



The "classical" (+) vertices have an off-diagonal struc- 
ture in Kedlysh space, while the quantum (-) vertices are 
diagonal 



v/2" 



-1! 



V2 



1. 



With this information, the conduction electron self en- 
ergy is given by 



£(*) = 



E«(t) Ejc(t) 
£ A (t) 



E+_(t) £++(<) 







-it) 



which in expanded form gives 



y ' 2N 



Tr[ n S(i)J(-t)] TrfnB^TiJH)] 
Tr[B(t)nJ(-t)] 



The Keldysh traces give the combinations 

Tr[nBj] =B r Jk + B k Ja 
Tr[ n B n J] = B R J A + B A J R + B K J K 

Tt[B Ti J]=B k J a +B r J k (79) 

Carrying out the Fourier transforms, the explicit form for 
the conduction self-energy is then 



£*(e) 

£ A ( e ) 
E*(e) 



B K (u)J A (uj-e)], 
B A (u;)J K (u-e)], 



= -it z 



K f duj 



B k (lo)J k (w-e) 



K i 



N J 4tt 

B a (u)J r (lo - e) + B r (lo) J a {u - e)] .(80) 



Repeating the same procedure for the slave fermions, we 
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obtain 



Similarly, 



+ B K (u)G A (u>-e)], 
n A ( £ ) = -*t 2 j^[B A >)G^-e) 
+ B a (uj)G k (ll> — e)] , 



II* (e) 



-it' 



du> 

4"7T 



B {w)G K (u - e) 



+ B a (lo)G r (lu - e) + B H (u;)G A (a; - e)] (81) 



Since the Schwinger boson self-energy vanishes in the 
large N limit, the Schwinger boson fields are unrenor- 
malized and take their bare values 



B R > A (v) = ' 
B K (y) = -2ni8{v - X)hi 



(82) 



where h B = 1 + 2n B and n B = n B (X) = 2SK/(NK) 
is the Bose distribution function. The Keldysh Green's 
function of the fermions can be related to their retarded 
and advanced Green functions via the relations 

J k {uj) = (J R (w)- J A {uo))h x {uj) 

= 2iIm[J R {u>)}(l-2f x (cj)), 

G k (uj) = (G R {u>) - G A (u))h*(u) 

= 2iIm[G R (u)](l-2fu(w)). (83) 



S^(e) = -A B (X)J R {X-e) 



^VJ-oo 7T 'uj-X + e-i5 
,K 



-ImJ R {uj) 



£*(e) = -t 2 -n B (X)J A {X-e) 



i 1 



K f°° dw 



N J~oc, i" 'uj-X + e + iS 
,K 



-ImJ H {uj) 



2n B (A)/ x (A-e)] 



S K (e) = -2it 2 -ImJ H (X-e)[n B (X)-f x (X-e) 



(87) 



We can also compute the electron distribution function 
from 



(w) 



E K h x h B - 1 



- * A h B - h x 



but this recovers exactly the result obtained in 
showing that detailed balance is satisfied. 



Appendix II: on the derivation of general current 
formula in case of two channels 



As an example we present the derivation of the Green's 



function G 



lip 



G ltp {t x M) = -i{T ck d x (ti)y + {t 2 )), (89) 
and the Dyson Eq is that depicted on Fig. 4 To avoid dou- 



With these simplifications, we can expand the self- 
energies in the following form 

IL A (e) = -t 2 n B (X)G R {X~e) 

- t 2 r -1 -M K H 

J^oo it ui — A + e — id 

IL R (e) = —t 2 n B (X)G A (X — e) 

/*H r— | ■ —IlTlG R (ld) 



cko 

w — A + e + 

n K (e) = -2ii 2 /mG fl (A-e)[n s (A)-/*(A-e) 

- 2n fl (A)/*(A-e)] (84) 

The slave fermion distribution function is given by 



11^ hyh B — 1 

lx = m - n A = hB-h*- 



(85) 



Rearranging this expression gives 

h x h B — 1 



h B — h x 



(86) 



G, 



< 



FIG. 9: The Dyson Eq. for Gi v Green's function (Gi v is a 
matrix in Keldysh space). 

ble counting we renormalize only the (i-electron Green's 
function, while the channel Green's function is left un- 
renormalized, so that the expression for G\ a is 

Gi v (ti,t 2 ) = / dr (D 11 (t 1 ,T)t 1 g^ v (T,t 2 ) 

JcK v 



where 



+D u {ti,T)t 2 g 9V (T,t 2 )) (90) 



Ai(*i,*a) = -iiTckdiit^dfih)). 



(91) 



Expression (|90[) can be rewritten in Larkin-Ovchinnikov 
space and the Keldysh part of the Green's function can 
be found. After Fourier transforming, we obtain the fol- 
lowing Keldysh component 

G^=ti(D R g^+D^g A 9 )+t 2 (D R g^+D^g A (p ) (92) 
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The Dyson Eq for G 2 ^ is derived in an analogous fash- 
ion 

G% = t2{DKg^+D«g^)Mx{D*g^+Dgg^) (93) 

These equations can now be expanded using the explicit 
expressions for the bare channel Green's functions given 
in section IV, B to obtain the general expressions (|56|) and 

Appendix III: estimation of the interference current 

term 



a) 



b) 




1 ^* ' 1 2 

FIG. 10: Diagrammatic contributions to the off-diagonal 
Green's function Di2(u>), describing inter-channel contribu- 
tions to the conduction electron propagator. 

Contributions to the current which are produced by 
channel interference are defined by the interchannel 
Green's functions D 12 and D 2 i- These Green's functions 
are proportional to the self-energies of G-functions. On 
Fig. 10 we show two types of contributions which arise 
in the first order of the perturbation theory. Diagrams 
of type a) contribute to the non-elastic scattering of d- 
electrons, and therefore constitute the non-Fermi liquid 
corrections to the conductance. Other diagrams can be 
summed in a simple Dyson equation 

Du = D n tiG^, 2 — D n tig^tiD 12 + D^tig^^D^, 

D12 = t x t 2 {l - tjD^^D^g^D^i) 



where, as usual are bare Green's functions, g^, g^ip- 
bare conduction electron Green's functions, and D\ 2 , D 22 
are the fully-renormalized d-electron Green's functions. 
Analogously we can write the Dyson Eq for Du 

D\\ — D n + tlD^g^Du + ti^D^g^^D^ 

= D^+tlDng^D^+hhDug^D ^ (95) 
We can write from here 

£>ii(l - t 2 i9^D n ) = 25^(1 + hhD 12 g^), (96) 

so that 

D u = (1 + hhD^g^D^l - tlg^D^)- 1 (97) 
and 

(1 + h^D^g^^Dn = D^l - tlg^D^)' 1 

= (l-t 2 1 D n g^)- 1 D° 11 (98) 

Now we can plug (|9"8|) into (|M|) and get 

D 12 = (1 + t 1 t 2 D 12 g^)- 1 t 1 t2D 11 g M , v D22 (99) 

which leads to the equation 

D 12 = t 1 t 2 D 11 g lPv D 2 2 - ht 2 D 12 g vij D 12 . (100) 

From (|100p we now can derive that retarded and ad- 
vanced off-diagonal terms vanish in this approximation. 
Namely, for example 

D* = t 1 t 3 D? 1 g$ v D& - h^DKg^D* (101) 

is correct and therefore D^ 2 = 0, because g5 = g^ — 0. 
Keldysh components of D\ 2 nevertheless do not vanish 

Df 2 = ht 2 (D?x4 v Dg+ + D^D*) + 

+ tMDKg^D* + D? 29 ^Df 2 + DKg^DfeW) 

Most of the terms in (| 1 02|) vanish and we are left with 

DF 2 =t 1 t 2 DKg« v D 2 A 2 . (103) 
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